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Abstract 

The size of the particle emission region in high-energy collisions can be deduced using the femto- 
scopic correlations of particle pairs at low relative momentum. Such correlations arise due to quan¬ 
tum statistics and Coulomb and strong final state interactions. In this paper, results are presented 
from femtoscopic analyses of K^K^, K^K^, pp, and pp correlations from Pb-Pb collisions 

at a/snn = 2.76 TeV by the ALICE experiment at the LHC. One-dimensional radii of the system 
are extracted from correlation functions in terms of the invariant momentum difference of the pair. 
The comparison of the measured radii with the predictions from a hydrokinetic model is discussed. 
The pion and kaon source radii display a monotonic decrease with increasing average pair transverse 
mass mj which is consistent with hydrodynamic model predictions for central collisions. The kaon 
and proton source sizes can be reasonably described by approximate mj-scaling. 
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1 Introduction 


Two-particle correlations at low relative momenta (commonly referred to as femtoscopy), which are 
sensitive to quantum statistics (in the case of identical particles) as well as strong and Coulomb final- 
state interactions (ESI), are used to extract the space-time characteristics of the particle-emitting sources 
created in heavy-ion collisions M] . The source radii extracted from these correlations describe the sys¬ 
tem at kinetic freeze-out, i.e. the last stage of particle interactions. Pion femtoscopy, which is the most 
common femtoscopic analysis, has shown signatures of hydrodynamic flow in heavy-ion collisions, man¬ 
ifesting as a decrease in the source radii with increasing transverse mass mx = \ -|- mf 0,11], where 


kj = Ipt.i -|-Pt,2|/2 is the average transverse momentum of the pair. This behavior can be interpreted 
as one of the signatures of the formation of deconfined quark matter in these collisions 10]. However, 
a necessary condition for collective behavior is for all particles created in the collision, not just pions, 
to experience hydrodynamic flow. Thus, femtoscopic studies with particles other than pions are also 
needed. It was shown that the hydrodynamic picture of nuclear collisions for the particular case of small 
transverse flow leads to the same mj behavior of the longitudinal radii (/?iong) for pions and kaons |0]. 
This common mx-scaling for n and K is an indication that the thermal freeze-out occurs simultaneously 
for 71 and K and that these two particle species are subject to the same Eorentz boost. Previous kaon 
femtoscopy studies carried out in Pb-Pb collisions at the SPS by the NA44, NA49, and CERES Collabo¬ 


rations Il8l4l0l] reported the decrease of Riong with mj as ~ as a consequence of the boost-invariant 


longitudinal flow. Subsequent studies carried out in Au-Au collisions at RHIC iliMl have shown the 
same level in the mx-dependencies for 7t and K radii, consistent with a common freeze-out hyper-surface. 
Eike the SPS analysis, no exact universal mx-scaling for the 3D radii was observed at RHIC. In the case 
of the one-dimensional correlation radius /?inv, only approximate scaling with nij is expected as an ad¬ 
ditional confirmation of hydrodynamic expansion 1^. In fact, source sizes as a function of for 
different particle types (tt, K, p...) follow the common curve with an accuracy of ~ 10%. 

The motivation for comparing femtoscopic analyses with different particle species is not limited to study¬ 
ing mx dependence. The kaon analyses also offer a cleaner signal compared to pions, as they are less 
affected by resonance decays, while the proton analysis provides a possibility for checking if baryons 
are included in the collective motion. Studying charged and neutral kaon correlations together provides 
a convenient experimental consistency check, since they require different detection techniques (charged 
tracks vs. decay vertex reconstruction) and call for different final-state interaction fitting parametriza- 
tions (Coulomb-dominated vs. strong interaction-dominated), yet they are predicted to exhibit the same 
femtoscopic parameters 11411 . In addition to the charged kaon analyses at the SPS and RHIC, neutral 
kaon correlations were studied in Au-Au collisions at RHIC s, and AEICE has performed analyses 
on both charged and neutral kaons in pp collisions EH- Recent pion femtoscopic results were ob¬ 
tained at RHIC 1 18] and the EHC Js, 3-21], and proton femtoscopy has also been previously studied at 
RHIC S- 


This paper presents the results of femtoscopic studies of K^K^, K^Kg, pp, and pp correlations 

from Pb-Pb collisions at .^^nn = 2.76 TeV by the AEICE experiment at the EHC. The femtoscopic radii 
and A parameters (the latter describe the decrease of the femtoscopic correlations due to e.g. long-lived 
resonances; see Sec. 13.1! and |4ll are extracted from one-dimensional correlation functions in terms of 
the invariant momentum difference for a range of collision centralities and mx values. A hydrokinetic 
model lll4l is used to compare the kaon experimental results with hydrodynamic predictions. 


The organization of the paper is as follows. In Sec. |2l we describe the data selection criteria. In Sec. [3j 
the details of the correlation functions and the fitting process are discussed. The results of the analysis 
are shown in Sec. HI and a summary is provided in Sec.[5j 


2 












One-dimensional femtoscopy in Pb-Pb collisions 


ALICE Collaboration 


2 Data analysis 


The dataset analyzed in this paper is from Pb-Pb collisions at -^^nn = 2.76 TeV at the LHC measured 
by the ALICE detector 12311 . About 8 million events from 2010 and about 40 million events from 2011 


were used (2010 data were analyzed in the pion and K® analyses only). Events were classified according 


to their centrality determined using the measured amplitudes in the VO detectors 12411 . Charged particle 
tracking is generally performed using the Time Projection Chamber (TPC) 12511 and the Inner Tracking 
System (ITS) 12311 . The ITS allows for high spatial resolution in determining the primary (collision) 
vertex. In the pion, charged kaon, and proton analyses, the determination of the momenta of the tracks 
was performed using tracks reconstructed with the TPC only and constrained to the primary vertex. 
Primary tracks were selected based on the distance of closest approach (DCA) to the primary vertex. 
Additional track selections based on the quality of the track reconstruction fit and the number of detected 
“hit” points in the TPC were used. Also, all primary pairs sharing more than 5% of TPC clusters were 
rejected. In the neutral kaon analysis, the secondary daughter tracks used global (TPC and ITS) track 
reconstruction and did not use any cuts based on track reconstruction quality or number of used or shared 
TPC clusters. The secondary vertex finder used fo locafe fhe neufral kaon decays employed fhe “on- 
Ihe-fly” reconsfrucfion mefhod 12611 . which recalculafes fhe daughter frack momenfa during fhe original 
fracking process under fhe assumption fhaf fhe fracks came from a decay vertex insfead of fhe primary 
vertex. 

Particle idenfificafion (PIP) for reconsfrucfed fracks was carried ouf using bofh fhe TPC and fhe Time- 
of-Elighf (TOP) defector 12711 in fhe pseudorapidify range |t]| < 0.8. Por TPC PIP, a parametrized 
Bethe-Bloch formula was used to calculate the specific energy loss (dE/dv) in fhe defecfor expected for 
a parficle wifh a given mass and momenfum. Por PIP wifh TOP, fhe parficle mass was used fo calculafe 
fhe expecfed fime-of-flighf as a funcfion of frack lengfh and momenfum. Por each PIP mefhod, a value 
N(y was assigned to each frack denoting fhe number of sfandard deviations befween fhe measured frack 
informafion and fhe calculations menfioned above. Pifferenl cuf values of Na were chosen based on 
defector performance for fhe various particle types and track momentum (see Table [T] for specific values 
used in each analysis) i28n . 


The analysis defails specific fo each parficle species used in fhis sfudy are discussed separafely below. 


2.1 Pion selection 

The main single-particle selecfion criferia used in fhe pion analysis are summarized in Table [T] Pion 
idenfificafion was performed using fhe TPC only. An overall purify of fhe pion candidate sample was 
esfimafed using TPC dE/dx disfribufions of fhe dafa and was found fo be above 95%. The main source 
of confaminafion comes from e^ in fhe region where fhe dE/dx curves for pions and elections infersecf. 

Pemfoscopic correlafion funclions of identical parficles are sensifive fo fhe fwo-frack reconsfrucfion ef¬ 
ficiency because fhe correlated particle pairs (i.e. fhose wifh small relative momenfum) generally have 
close tiajecfories. The main fwo-frack issues are spliffing (fwo fracks reconstructed from one particle) 
and merging (one track reconstructed from two particles), which are generally avoided using a track sepa¬ 
ration cut. Por pions, pairs were required to have a separation of |At] | > 0.016 or y^Aq^ > 0.045 

measured at the radial distance 1.2 m. Here, rj is the pseudorapidity, and (p* is the azimuthal coordinate 
taking into account track bending due to the magnetic field. 


2.2 Charged kaon selection 

The main single-parficle selecfion criferia used in fhe charged kaon analysis are listed in Table [T] K* 
idenfificafion was performed using fhe TPC (for all momenfa) and TOP (for p > 0.5 GeV/c) defecfors. 
Pigure [T] (a) shows fhe momenfum dependence of fhe single kaon purify, defined as fhe fracfion of ac- 
cepfed kaon fracks fhaf correspond fo tiue kaon parficles. The purify values were obfained from TPC 
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Pion selection 


Transverse momentum pj 

0.14 < pt <2.0 GeV/c 

\n\ 

<0.8 

Transverse DCA to primary vertex 

< 0.2 cm 

Eongitudinal DCA to primary vertex 

< 0.15 cm 

A^ct.TPC 

<3 


Charged kaon selection 


Pt 

0.15 <pt < 1.5 GeV/c 

1^71 

<0.8 

Transverse DCA to primary vertex 

< 2.4 cm 

Eongitudinal DCA to primary vertex 

< 3.0 cm 

TVff.TPC (for P < 0-5 GeV/c) 

<2 

A^ct.tpc (for p > 0.5 GeV/c) 

<3 

TVct.tof (for 0.5 < p < 0.8 GeV/c) 

<2 

A^ff.TOF (for 0.8 < p < 1.0 GeV/c) 

< 1.5 

A^ct.tof (for 1.0 < p < 1.5 GeV/c) 

< 1.0 


Neutral kaon selection 


1^1 

< 0.8 

Daughter-daughter DCA 3 D 

< 0.3 cm 

DCA 3 D to primary vertex 

< 0.3 cm 

Decay length 

< 30 cm 

Cosine of pointing angle 

>0.99 

Invariant mass 

0.480 < < 0.515 GeV/c^ 

Daughter px 

>0.15 GeV/c 

Daughter Tj 

<0.8 

Daughter DCA 3 D to primary vertex 

> 0.4 cm 

Daughter Na jpc 

<3 

Daughter Na, tof (for p > 0.8 GeV/c) 

<3 


Proton selection 


Pt 

0.7 < Pt < 4.0 GeV/c 

1^1 

<0.8 

Transverse DCA to primary vertex 

< 2.4 cm 

Eongitudinal DCA to primary vertex 

< 3.2 cm 

Nff jpc (for p < 0.8 GeV/c) 

<3 

\/^ct,tpc +^ct,tof (foJ" P > 0-^ GeV/c) 

<3 


Table 1: Single particle selection criteria. 


AEldx distributions of the data and by studying HIJING simulations using GEANT3 to model 
particle transport through the detector. Eike the pions, the dominant contamination for charged kaons in 
the momentum region 0.4 < p < 0.5 GeV/c comes from e^. The pair purity is calculated as a product 
of two single-particle purities, where the momenta are taken from the experimental distribution. The 
pair purity as a function of kj for three different centralities is shown in EigiHb). 

Regarding two-track selection criteria, charged kaon pairs were required to have a separation of | At] | > 0.02 
or |A(p*| > 0.017 measured at the radial distance 1.6 m. 
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Fig. 1: (Color online) Single purity (a) and pair purity (b) for different centralities. In (b) the values for 
different centrality intervals are slightly offset for clarity. 


x10® 



Fig. 2: (Color online) Invariant mass distribution of tt+tt pairs showing the peak for two centrality 
intervals. The 45-50% centrality is scaled so that both distributions have the same integral in the range 
0.480 < < 0.515 GeV/c^. 


2.3 Neutral kaon selection 

The decay channel Kg —)• 71+was used for the identification of neutral kaons. The single-particle cuts 
for parents (Kg) and daughters (tt^) used in the decay-vertex reconstruction are shown in Table [T] PID 
for the pion daughters was performed using both TPC (for all momenta) and TOE (for p > 0.8 GeV/c). 
Eigure|2] shows an example of the tt+tt^ invariant mass distribution where the Kg peak is seen. The cuts 
used in this analysis were chosen to balance statistics and signal purity. The neutral kaon purity (defined 
as Sig./[Sig.-i-Bkg.] for 0.480 < rrij^+j^- < 0.515 GeV/c^) was found to be greater than 0.95. 
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Two main two-particle cuts were used in the neutral kaon analysis. To resolve two-track inefficiencies 
associated with the daughter tracks, such as the splitting or merging of tracks discussed above, a sepa¬ 
ration cut was employed in the following way. For each kaon pair, the spatial separation between the 
same-sign pion daughters was tabulated at several points throughout the TPC (every 20 cm radially from 
85 cm to 245 cm) and averaged. If the average separation of either pair of tracks was below 5 cm, the 
kaon pair was not used. Another cut was used to prevent two reconstructed kaons from using the same 
daughter track. If two kaons shared a daughter track, one of them was cut using a procedure which com¬ 
pared the two Kg candidates and kept the candidate whose reconstructed parameters best matched those 
expected of a true Kg particle in two of three categories (smaller Kg DCA to primary vertex, smaller 
daughter-daughter DCA, and Kg mass closer to the PDG value Olll l. This procedure was shown, using 
HIJING-I-GEANT3 simulations, to have a success rate of about 95% in selecting a true Kg particle over 
a fake one. More details about KgKg analysis can be found in Ref. 


2.4 Proton selection 

The single-particle cuts used in the proton analysis are summarized in Table [T] The proton analysis used 
tracks with 0.7 < px < 4.0 GeV/c. The lower px cut is used to suppress protons coming from weak 
decays and interactions with the detector material. Particle identification for p and p was performed 
using both TPC (for all momenta) and TOE (for p > 0.8 GeV/c). The proton purity was estimated 
using HIJING-I-GEANT3 simulations and was found to be greater than 95%. The used DCA criteria 
do not fully discriminate between primary protons and protons from weak decays. This may lead to 
a significant contamination from protons from lambda particles. The effect of this contamination is 
discussed in Sec. 13.41 

Regarding two-track selection criteria, pairs were required to have a separation of |At] | > 0.01 or |A(p*| > 
0.045 measured at the radial distance 1.2 m. 


3 Construction of the correlation functions and fitting procedures 


The experimental two-particle correlation function is defined as C(q) = A(q)/R(q), where A(q) is the 
measured distribution of same-event pair momentum difference, q = Pi — P 2 ^ and R(q) is the reference 
distribution of pairs from mixed events. The pairs in the denominator distribution B(q) are constructed by 
taking a particle from one event and pairing it with a particle from another event with a similar centrality 
and primary vertex position along the beam direction. Each event is mixed with five (ten) others for the 
Kg (tt^, K^, p) analysis. The available statistics of proton pairs with low q (< 0.2 GeV/c) allowed us to 
perform the analyses only for the one-dimensional correlation function C{q), where ^ = |q| in the Pair 
Rest Frame (PRF). In the case of pions and kaons, the statistics were high enough for three-dimensional 
studies, but these are beyond the scope of this paper; here, only the one-dimensional analysis is presented 
in order to compare results with heavier particles. The numerator and denominator are normalized such 
that C(^) —>• 1 as <7 —)■ oo. Pair cuts have been applied in exactly the same way for the same-event (signal) 
and mixed-event (background) pairs. 

All correlation functions have been corrected for momentum resolution effects. The correction factors 
were determined using HIJING events to build simulated correlation functions using theoretical corre¬ 
lation functions as weights. The ratio of the correlation functions using HIJING generated momenta to 
those using HIJING-I-GEANT3 reconstructed momenta forms the correction factor. 


For the analyses presented in this paper, the theoretical femtoscopic correlation function is defined as the 
square of the two-particle wavefunction averaged over the relative distance r* of the emitters in the PRF. 
This is performed using the Koonin-Pratt equation 133113411 


C{q)= j S{r*)\^>{q,r*)\^A\\ 


( 1 ) 
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For the one-dimensional analysis, we assume a spherically symmetric Gaussian distribution of the parti¬ 
cle emitter spatial separation r* in the PRF with size Rim 0511 . 


S(r*)~exp(-r* /4R^m)- 


( 2 ) 


The two-particle wavefunction is (anti-)symmetrized for identical bosons (fermions) and may include 
terms incorporating Coulomb or strong final-state interactions, depending on the type of particles being 
studied. 

The methods used in constructing and fitting the various correlation functions are discussed separately 
below. 


3.1 Pions 

Pion correlation functions were fitted using the Bowler-Sinyukov formula 

C{q) =N[\-X+XK{q) (l +exp , 


(3) 


where N is the normalization factor. The ^^arameter (also used in the other analyses) can be affected by 
long-lived resonances, coherent sources |13, ^ and non-Gaussian features of the particle-emission 
distribution. K{q) is a symmetrized A'-factor calculated according to Ref. iflil ls^ as 


K{q) = C(QS + Coulomb)/C(QS), 


(4) 


where C(Q^and C(QS -I- Coulomb) are the theoretical correlation functions calculated with THERMI- 
NATOR 2 idOn using the quantum statistics (“QS”) and “QS-i-Coulomb” weights (i.e. squared wavefunc¬ 
tion), respectively idlH . The effect of the strong interaction is neglected here, since for like-sign pions, 
the contribution is small for the expected source sizes idlH . Figure [3] shows an example 71+71+ correla¬ 
tion function with the corresponding line of best fit. More details about the pion analysis may be found 
in llH. 


3.2 Charged kaons 

Figure |4] shows an example K+K+ correlation function with the corresponding line of best fit. A purity 
correction was applied to the correlation function according to 


Ccorrected — (Craw 1 TP)/P, 


(5) 


where the pair purity P is taken from Fig. [T] K+K+ correlation functions were fitted using the Bowler- 
Sinyukov formula of Eq. [3l the procedure is essentially the same as for pions. There are no available 
experimental data for K+K+ strong FSI. The influence of the strong interaction to the correlation function 
was estimated with the s-wave scattering length calculated within the fully-dynamical lattice QCD 14311 . 
The systematic uncertainty assigned to this effect was determined to be 4%. 

3.3 Neutral kaons 

Figure [5] shows an example K^Kg correlation function with the corresponding line of best fit. K^Kg 
correlation functions were fitted with a parametrization which includes Bose-Einstein statistics as well 
as strong final-state interactions (FSI) ilSL 13511 . 


C{q) =[l-X+XC'{q)]{a + bq), (6) 

where 

C\q) = 1 + + CstrongFSl(^,^), (V) 
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Fig. 3: (Color online) Example correlation function with fit for 7r+7r+ for centrality 5-10% and (k^) = 0.35 GeV/c. 
Statistical uncertainties are shown as thin lines. 



Fig. 4: (Color online) Example correlation function with ht for for centrality 0-10% and (A:t) = 0.35 GeV/c. 

Systematic uncertainties (boxes) are shown; statistical uncertainties are within the data markers. The main sources 
of systematic uncertainty are the momentum resolution correction and PID selection. 
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Fig. 5: (Color online) Example correlation function with fit for K^Kg for centrality 0-10% and (A:t) = 0.48 GeV/c. 
Statistical (thin lines) and systematic (boxes) uncertainties are shown. The main source of systematic uncertainty 
is the variation of single-particle cuts. 



Fig. 6: (Color online) Example correlation function with ht for pp for centrality 0-10% and (kj) = 1.0 GeV/c. 
Statistical (thin lines) and systematic (boxes) uncertainties are shown. The main source of systematic uncertainty 
is the variation of two-track cuts. 
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CstrongFSI ) — 


1 


f{q) 


R 


+- {qR) - ■^F 2 {qR) 


^/nR 


R 


and 


Fi(z)= / dx-; F2{z) = - 

Jo z z 


( 8 ) 


(9) 


f{q) is the s-wave scattering amplitude for the system; we neglect the scattering for and 

due to small scattering lengths ss 0.1 fm 115]. The factor of 1/2 in Eq. [^is due to the fact that 
half of the KgKg pairs come from The strong ESI have a significant effect on the con¬ 
tribution to the correlation function due to the near-threshold resonances, fo(980) and ao(980). 

Eor the scattering amplitude, only s-wave contributions were taken into account; the higher-order cor¬ 
rections were small and therefore neglected 14411 . The scattering amplitude f{q) is calculated using a 
a two-channel parametrization which accounts for the elastic transition —)■ and the inelas¬ 
tic transition K+K^ —)• (See Ref. ilSll for more detailed expressions describing the fit function). 

Equation 12 also includes an additional factor to account for non-femtoscopic background correlations at 
large q, with a and b being free parameters in the fit. 


3.4 Protons 


Eigure [2 shows an example pp correlation function with the corresponding line of best fit. The femto- 
scopic correlations of pp and pp pairs are due to a combination of Eermi-Dirac statistics. Coulomb and 
strong ESI. A distinct maximum is seen at ^ « 40 MeV/c iSSH ; this enhancement is due to the strong 
interaction, as both quantum statistics and Coulomb interaction present a negative correlation. Due to 
the fact that feed-down from weak decays cannot be neglected in high-energy heavy-ion collisions, the 
effects of residual correlations related to the pA system are taken into account. The proton daughter of 
a A decay has similar momentum to the A itself and may survive the experimental selection for primary 
protons. Thus, it may contribute to the measured correlations by forming a pair with a primary proton. 
As can be seen in Eig. [2 attempting to fit the measured correlation functions with the theoretical pp 
(pp) functions alone was unsuccessful due to the additional positive correlation observed in the range 
60 < ^ < 160 MeV/c. Thus, a method of simultaneous fitting of pp (pp) and pA (pA) correlations was 
applied. Contributions from heavier baryon-baryon pairs are not taken into account since the original 
correlation between the parent particles is not known due to unknown interaction parameters, for exam¬ 
ple for the AA pair. Moreover such residual correlations are more smeared compared with pA because 
of larger decay momentum. In addition, the fraction of baryons heavier than A decaying to protons is 
smaller than the fraction of A’s. Einally, comparing with baryon-antibaryon pairs analysed in ll45n . the 
width of the correlation for baryon-baryon pairs is much smaller, and therefore the effect is much more 
smeared due to decay kinematics. 


The experimental correlation function of pp and pp systems were fitted with 1451] 


Cmeas(^pp) — 1 T ^p( 6 ^pp(^pp 5 ^) f) T ^A(CpA(^pp 5 /^) 1)) (10) 


where App is the fraction of correlated pp pairs where both particles are primary, and ApA is the fraction of 
correlated pp pairs where one particle is primary and the other is a daughter of A decay. The theoretical 
proton-proton correlation function was calculated as 


C'(^pp) 


1 


3 


4 

fS(r*) 

+ 4 

fS(r*) 


( 11 ) 


This formulation takes into account the necessary (anti-)symmetrization of the wavefunction for a pp 
pair in the singlet (triplet) spin state with a corresponding weight of 1/4 (3/4). The pp pair wavefunction 
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may be written as 1441] 






G(P,V) 


( 12 ) 


where r* is the spatial separation of particle emitters at generally different emission moments in the 
PRF, Sc = argr(l +iT]) is the Coulomb s-wave phase shift, Ac(ri) = 27iri{e^^^ “ 1) ^ is the Gamow 
factor (also referred to as Coulomb penetration factor), T] = (^a^pp)^', a = is the two- 

particle Bohr radius taking into account the sign of the interaction (a = 57.6 fm for pp pair), F is the 
confluent hypergeometric function, ^ -|-cos d*), 0* is the angle between Qpp and r*, G is the 

combination of the regular and singular s-wave Coulomb functions, and p = ^^ppr*. The amplitude of 
the low-energy s-wave elastic scattering due to the short range interaction /c(^pp) may be expressed as 


/c(^pp) 


1 ‘5io<?pp 

/o 8 


1 2 
-i^ppAc(Tj)--/j(?7) 


n -1 


(13) 


where /o is the scattering length, do is the effective radius of the interaction, h{ri) = [t/r(iT]) -|- t/r(—ip) — 
ln(p^)]/2, and xj/ is the digamma function. For the pp system in the singlet (triplet) state, fo and do are 
7.77 fm (-5.4 fm) and 2.77 fm (1.7 fm). 

For the feed-down term, the theoretical pA correlation function for a given Rpj^ transformed into the pp 
momentum space is obtained from the Fednicky-Fyuboshitz model 13511 and calculated as 


CpA(^pp>RpA) — ^C'pa(Ppa>7?pa)F(<7pp,<7pa)/^ T(ppp,Ppa), 


^pA 


^pA 


(14) 


where CpA(ppA;RpA) = 1 +CstrongFSi(^pA;^pA )5 and T{qp„^qpA) are the transformation factors related to 
A decay kinematics, calculated with THERMINATOR 2 14011 . Here, a spin-dependent version of Eq.[8]is 
used 


CstrongFSl(^)-^) — 

S 


fid) 


R 


1 - 


% 


^ m^jq) 

2y/nRJ ' ^/KR 


+ 


FiiqR)- 


^fjq) 

R 


FiiqR) 


(15) 


where f{q) is the spin-dependent scattering amplitude, is the fraction of pairs in each total spin state 
S, and Jq is the effective radius of the interaction. It is assumed that the radii of pp and pA sources are 
equal. Therefore, there are three free fit parameters in Eq. [T0| App, ApA , and R. Theoretical pp and 
pA correlation functions were calculated using several values of the free parameters, and the fit function 
(for the set of parameters given during each fit iteration) was formed by a quadratic interpolation of the 
calculated correlation functions. 


3.5 Systematic uncertainties 

The effects of various sources of systematic uncertainty on the extracted fit parameters were studied as 
functions of centrality and kx- Table[2]shows the minimum and maximum uncertainties from each source. 
The values of the total uncertainty are not necessarily equal to the sum of the individual uncertainties, as 
the latter can come from different centrality or kj bins. All four analyses studied the effects of changing 
the selection criteria for the events, particles and pairs used (variation of cut values up to ±25%) and 
varying the range of q values over which the fit is performed (variation of q limits up to ±25%). Un¬ 
certainties associated with momentum resolution corrections are included in the K, K^, and p analyses; 
the Kg analysis also studied this and found the uncertainties to be negligible. The K^, Kg, and p anal¬ 
yses encountered uncertainties associated with the non-flat background seen at large-^ for high-kj pairs 
in peripheral collisions (estimated by using different parametrizations (linear or polynomial) to fit the 
large-g region). Strong FSI uncertainties affect both kaon analyses. For Kg, the strong FSI uncertainty 
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71^ 


K± 



3 


^inv 

A 

^inv 

A 

^inv 

A 

^inv 

A 

Evenf/parficle/pair selection 

2-13 

6 

3-5 

3-5 

1-4 

2-10 

2-27 

12-58 

Non-Hal background 

- 

- 

0.2-5 

0.2-5 

0-5 

0-4 

2-3 

1-9 

Eil range 

10 

33 

1-5 

1-5 

0-4 

0-3 

3-26 

3-57 

Sfrong ESI 

- 

- 

4 

4 

1-2 

5-10 

n/a 

n/a 

Coulomb function 

3 

3 

2 

4 

n/a 

n/a 

n/a 

n/a 

PID and purify 

- 

- 

- 

5-28 

- 

- 

4-18 

13-41 

Momenlum resolution 

2 

3 

3-5 

5-10 

- 

- 

3 

1-8 

Eixing App 

n/a 

n/a 

n/a 

n/a 

n/a 

n/a 

1-29 

n/a 

Rpp/RpA ratio 

n/a 

n/a 

n/a 

n/a 

n/a 

n/a 

1-13 

20-52 

Total (quad, sum) 

11-21 

34 

6-9 

10-32 

2-7 

7-15 

10-40 

30-80 


Table 2: Minimal and maximal uncertainty values for various sources of systematic uncertainty (in percent). The 
A for the proton analysis refers to the sum of App and ApA. Please note that each value is the minimum (maximum) 
uncertainty from a specihc source, but each can be from a different centrality or kj bin. Thus, the minimum (max¬ 
imum) total uncertainties are greater (smaller) than (or equal to) the sum of the minimum (maximum) individual 
uncertainties, ’n/a’ denotes that the given descriptor of the systematic uncertainty is not applicable for the specific 
pair type, and means that the contribution from the given source is negligible. 


comes from the fact that several sets of fo(980) and ao(980) parameters are available l3-49]; each set 
is used to fit the data, the results are averaged, and the maximum difference was taken as the systematic 
error. The n and analyses have uncertainties associated with the choice of the Coulomb function 
used in the fitting procedure. The analysis had additional uncertainties due to the misidentification 
of particles and the associated purity correction. The p analysis also had uncertainties associated with 
the uncertainty in the ratio and attempts to fix fhe App paramefer using fhe single-parficle purify. 

All of fhe analyses were performed separafely for fhe fwo differenl signs of fhe ALICE dipole magnetic 
field, buf fhe resulfing systematic uncerfainfy was found fo be negligible in all cases. 


Systematic uncerfainfies on correlafion functions (Eigs. IH |5] and O were derived from fhe variation of 
single- and fwo-parficle cufs. 


4 Results 


Eigures |7] and [8] presenf fhe exfracfed fif parameters from K*K^, K^K^, and pp correlations 

for several infervals of cenfralify and fransverse mass. Bofh sfalislical and systematic uncerfainfies are 
shown. The qualify of fhe fifs used fo exfracf fhe shown paramefers can be assessed using fhe x'^/NDF 
values, which are in fhe ranges of 1.2-5.0, 0.8-3.5, 0.6-1.5, and 0.8-3.2 for fhe pion, charged kaon, neufral 
kaon, and proton analyses, respectively. 


Eigure |7] shows fhe exfracfed A paramefers vs. mj for several cenfralifies. The proton A is fhe sum of 
App and ApA from Eq. [TO] The values for all species measured lie mosfly in fhe range 0.3-0.7 and show 
no significanl cenfralify dependence. The values of A are less fhan unify due fo long-lived resonances 
which dilufe fhe correlafion functions and also lead fo non-Gaussian shapes of fhe correlafion functions, 
especially in fhe one-dimensional case lEoll . Resulfs for kaons and protons are consisfenf wifh each ofher 
af similar ifiT. Values of A for pions are lower fhan for kaons due to fhe sfronger influence of resonances; 
an additional cause could be a partial coherence of pions il9l] . 


Eigure [8] shows fhe exfracfed 7?inv paramefers vs. nij for several cenfralifies. Eor overlapping nij, fhe 
radius paramefers are mosfly consisfenf wifh each ofher wifhin uncerfainfies, fhough fhe pion radii are 
generally larger fhan fhe kaon radii. The Kg radii are slighfly higher fhan radii for cenfral collisions, 
buf fhe difference is less fhan fhe sysfemafic uncerfainfies. The radius paramefers show increasing size 


12 

























One-dimensional femtoscopy in Pb-Pb collisions 


ALICE Collaboration 


ALICE Pb-Pb = 2.76 TeV 


0-10% 10-30% 30-50% 

7t 

+ 

* 

☆ 

K 

-K- ■ 

□ 

★ 

K 

sKs ▼ 

V 

0 

pp A 

A 

♦ 

PP • 

O 




O 




t 1 



f 



_^_1_^^_ 

_^_ 

_^_ 

_1_ 


0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

(m^) (GeV/c^) 


Fig. 7: (Color online) A parameters (App-|-ApA in case of (anti)proton pairs) vs. nij for the three centralities 
considered for K^K^, K^K^, pp, and pp. Statistical (thin lines) and systematic (boxes) uncertainties are 

shown. The nij values for different centrality intervals are slightly offset for clarity. 



Fig. 8: (Color online) Ri^ parameters vs. mj for the three centralities considered for K^K^, KgKg, pp, 

and pp. Statistical (thin lines) and systematic (boxes) uncertainties are shown. 
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Fig. 9: (Color online) Comparison of the HKM model (see text) with measured kaon X (a) and Ri^v (b) parameters 
for 0-5% centrality. Statistical (thin lines) and systematic (boxes) uncertainties are shown. 


ALICE Pb-Pb = 2.76 TeV, 0-5% 

+ 7t±7l± 

■ KtK* 


T K^K^^HKM K*K* 
• pp oHKMpp 


_ (b) t 

J_^_I_I_^_I_^_I_L 


1 -4 ^ ALICE Pb-Pb t = 2.76 TeV, 0-5% 


1 . 2 - 
1 - 

0 . 8 - 

■ 

0 . 6 - 

0.4 -+ + 

0.2-(3) 


NN 
+71*71* 

-K"K" 

^HKM K-K- 



+ + 


+ + 


with increasing centrality as would be expected from a simple geometric picture of the collisions. They 
also show a decreasing size with increasing mj as would be expected in the presence of collective radial 
flow 10]. Both of these dependences can be seen in previous femtoscopic measurements i,i and 

also reinforce the interpretation that collective flow is present in these collisions for pions, kaons (neutral 
and charged), and protons alike. Deviations from exact mx-scaling of 7?inv can be explained as a conse¬ 
quence of the increase of the Lorentz factor with decreasing particle mass. In a hydrodynamic model 
scaling is observed for the three-dimensional radii measured in the Longitudinally Co-Moving System 
(LCMS). The transformation from LCMS to PRF involves a boost along the outward direction only, 
where the boost value is proportional to the transverse velocity of the pair and inversely proportional to 
the particle mass (for similar mx). Thus, a smaller mass leads to an increase in the boosted 7?out and, sub¬ 
sequently, /?inv in the PRF. Indeed, we observe such an effect in the data, as pion radii are systematically 
higher than kaon radii at the same mx. 


A comparison of a hydrodynamic flow -i- kinetics model, HKM lll4l1 . with the measured and X 
parameters for 0-5% centrality is shown in Fig. |9l The HKM values in Fig. |9]are specifically from K*K^, 
but the predictions for KgK® and K^K^ are consistent with each other. For Rmy, the charged kaon data 
show very good agreement with the predictions. The experimental data for the neutral kaons are again 
slightly higher than for the charged kaons, but this difference is still within systematic uncertainties. For 
A, both sets of kaon data match the decreasing trend with increasing kj exhibited by the HKM points, 
but the model slightly overpredicts the data. It is shown in il4ll that the most important resonances 
for KK pairs, K*(890) and 0(1020), do not significantly influence the X parameter (due to their low 
contribution), and the discrepancy between the model and experimental data can be explained by the 
lower experimental kaon purity and deviations of the experimental correlation function shape from a 
Gaussian distribution. For protons, the HKM prediction is compatible with the data. HKM calculations 
for one-dimensional pion radii are currently not available, but three-dimensional radii were reasonably 
reproduced by this model 15 ill . 


5 Summary 

Results from femtoscopic studies of K^K^, K^Kg, pp, and pp correlations from Pb-Pb collisions 

at .^5 nn = 2.76 TeV with ALICE at the LHC have been presented. The femtoscopic radii and X pa¬ 
rameters were extracted from one-dimensional correlation functions in terms of the invariant momentum 
difference. It was found that the emission source sizes of kaons and protons measured in these collisions 
exhibit transverse mass scaling within uncertainties, which is consistent with hydrodynamic model pre¬ 
dictions assuming collective flow. The deviation from the scaling for the pions can be explained as a 


14 








One-dimensional femtoscopy in Pb-Pb collisions 


ALICE Collaboration 


consequence of the increase of the Lorentz factor with decreasing particle mass during the transforma¬ 
tion from LCMS to PRF systems 115On . The extracted A parameters are found to be less than unity, as 
is expected due to long-lived resonances and non-Gaussian correlation functions. The predictions of the 
hydrokinetic model (HKM) for the one-dimensional femtoscopic radii for charged and neutral kaons and 
protons coincide well with the observations. 
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